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ABSTRACT 

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive 
method for observing the build-up of cosmological structure. In this paper we develop and 
test a new statistic based on anisotropies in the measured galaxy power spectrum, which is 
independent of galaxy bias and matches the matter power spectrum shape on large scales. The 
amplitude provides a constraint on the derivative of the linear growth rate through /<x g (mass). 
This demonstrates that spectroscopic galaxy surveys offer many of the same advantages as 
weak lensing surveys, in that they both use galaxies as test particles to probe all matter in 
the Universe. They are complementary as redshift-space distortions probe non-relativistic ve- 
locities and therefore the temporal metric perturbations, while weak lensing tests the sum 
of the temporal and spatial metric perturbations. The degree to which our estimator can be 
pushed into the non-linear regime is considered and we show that a simple Gaussian damping 
model, similar to that previously used to model the behaviour of the power spectrum on very 
small scales, can also model the quasi-linear behaviour of our estimator. This enhances the 
information that can be extracted from surveys for ACDM models. 

Key words: cosmology: large-scale structure 



1 INTRODUCTION 



It is now widely accepted that the large-scale structure we see in 
the distribution of galaxies arises through a process of gravitational 
instability, which amplifies primordial fluctuations laid down in the 
very early Universe. The rate at which structure grows from these 
small perturbations offers a key discriminant between cosmological 
models. For instance, dark energy models in which general relativ- 
ity is unmodified predict different Large-Scale Structure formation 
compared with Modifie d Gravity models with the same background 
expansion (e.g. iDvali etaLlbOOOl: [Carroll et al]|2004l : lBranj|2005l ; 
iNesseris & Perivolaropoulosll2008l) . 

Structure growth is driven by the motion of matter (and in- 
hibited by the cosmological expansion). Galaxies are expected to 
act as test particles within this matter flow, so the motion of galax- 
ies carries an imprint of the rate of growth of large-scale structure. 
Because of this, many previous analyses have shown that obser- 
vations of these galax y peculiar velocities can distinguish between 
classes of models (e.g Jjain & Zhanj|2007l ; ISong & Kovama!l2008l ; 
ISong&Percivalll2008li . A key technique to statistically measure 
the growth of the v elocity field u ses redshift-space distortions seen 
in galaxy surveys jKaiseill 19871) . Galaxy maps, produced by esti- 
mating distances from redshifts obtained in spectroscopic galaxy 
surveys, reveal an anisotropic galaxy distribution. The anisotropies 
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arise because galaxy recession velocities, from which distances are 
inferred, include components from both the Hubble flow and pecu- 
liar velocities from the comoving motions of galaxies. These dis- 
tortions encode information about the build-up of structure. 

Many previous surveys have been analysed to measure /3 ~ 
Q.®f/b, where b is the deterministic, local, linear bias of the galax- 
ies. The latest generation of large surveys have provided ever 
tighter constraints. An alyses using the 2-d egree Field Galaxy Red- 
shift Survey (2dFGRS; Colless et al. 2003) have mea sured redshift- 



space distortions in both th e correlation function ( Peacock et al 



l200ll ; lHawkins et alj 120031) and power spectrum (Per cival et al 
| 2004|) . Using the Sloan Digital Sky Survey (SDSS; lYork et al 
2000), redshift-space d istortions have also been measured in 
the correlation function JZehavi et al]|2005l ; lOkumura et alj|20"08l ; 
ICabre & Gaztafiag j |2008), and using an Ei genmode decomposi- 
tion t o separate real and redshift-space effects iTegmar k et alj2004l 
1 20061) . These studies were recently extended to 7. — 1 Jguzzo etal 



2007) using the VIMOS-VLT Deep Survey (VVDS; Le Fevre et al 



2005; Garilli et alj|2008h and the 2SLAQ s urvey (Ida A ngela et al. 
2008). In addition to measuring j3 at z = 0.8. lLe Fevre et alj b005T) 
emphasised the importance of using large-scale peculiar velocities 
for constraining models of cosmic acceleration. 

On linear scales the theor y behind the observed re dshift-space 
distortions is well developed (Kaiser 1987; Hamilton 1998j), and 
previous analyses have made use of this theory to measure ji. On 
quasi-linear and non-linear scales we are instead reduced to mak- 
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ing approximations, or using fitting formulae based on numerical 
simulations. The standard approach is to use a 'streaming' model, 
where linear theory is spliced together with an approximation for 
random motion of particle s in co llaps ed obje c ts (see section O. 
iTinker. Weinberg & Zhengl d2006h and iTinkeij (2007), which are 
developments of work in lHatton&ColeUl999l) . discuss a number 
of possible improvements to the streaming model based on fits to 
numerical simulation results. These approaches aim to allow us to 
extract information on p from scales where the power spectrum 
does not match its linear form. One concern is that the theoretical 
dependence on ft might change in the quasi-linear regime, leading 
to complicated dependencies and that fits to simulations will only 
be correct in a subset of cosmologies and galaxy formation models 
similar to those from which the fits were derivecfl 

A better approach would be to analyse the physics behind 
redshift-space distortions, and to try to use this to devise the best 
estimator with which to extract cosmological information. Because 
galaxy velocities only depend on the distribution of matter, we can 
devise an estimator P(k) that is not affected by galaxy bias: instead 
it measures the large-scale shape of the matter power spectrum. 
Because the estimator is based on the velocity power spectrum, the 
large-scale amplitude is f 2 times that of matter fluctuations, where 
/ is the logarithmic derivative of the linear growth rate with re- 
spect to the scale factor, / = dlnD/dlna. From this estimator, we 
can measure /0"8, maS s, which is proportional to dD/d In a, and pro- 
vides a good discriminant between modified gravity and dark en- 
ergy models dSong & Percivall2008l) . P(k) is expected to match the 
linear matter power spectrum shape on large scales, so small-scale 
differences will help us to determine the scales on which the the- 
ory is breaking down. For ACDM models, we show that a simple 
Gaussian smoothing of the redshift-space power along the line-of- 
sight is able to match most of this break-down. This places redshift- 
space distortion measurements on the same footing as weak lensing 
measurements in the sense that they both allow us to test the mat- 
ter distribution directly. They provide complementary information, 
as non-relativistic velocity measurements only depend on tempo- 
ral metric perturbations, while weak lensing tests the sum of the 
temporal and spatial metric perturbations. 

The layout of our paper is as follows. We first review the 
theory behind redshift-space distortions (Section O in the linear, 
quasi-linear and non-linear regimes. In Section [3] we introduce a 
new estimator, based on the monopole and quadrupole from a Leg- 
endre decomposition of the redshift-space power spectrum, which 
is designed to recover a power spectrum given by f 2 times the mat- 
ter power spectrum on large scales. Section [4] introduces the N- 
body simulation used to test this theory and our new estimator. We 
use spherically averaged power spectra where we include distor- 
tions along multiple axes to analyse this simulation, an approach 
described in Section [5] Our analytic theory is compared with re- 
sults from the numerical simulation in Section [6] The paper ends 
with a discussion of our results. 



where u z (x) is the line-of-sight component of the galaxy veloc- 
ity (assumed non-relativistic) in units of the Hubble velocity, and 
we have taken the line-of-sight to be the z-axis. We shall adopt 
the "plane-parallel" approximation, so this direction is fixed for all 
galaxies. 

The galaxy overdensity field in redshift-space can be obtained 
by imposing mass conservation, (1 + S g )d 3 s = (1 + 6 g )d 3 r, and the 
exact Jacobian for the real-space to redshift-space transformation 



d 3 r 



2 ' du. 



(2) 



In the limit where we are looking at scales much smaller than the 
mean distance to the pair, ujz is small and it is only the second term 
that is important dKaiseJll98l but see IPapai & Sz apudi 200j|, 



l+<5;; = (l+<5 g )(l + 



du ; 
~dz 



(3) 



If we assume an irrotational velocity field we can write u- = 
<9/<9zV~ 2 6>, where 6 = V • u, and V~ 2 is the inverse Laplacian oper- 
ator. In Fourier space, (<9/<9z) 2 V~ 2 = (k,/k) 2 = //, where ji is the 
cosine of the line-of-sight angle, so we have that 



Sl(k) = 6 g (k) - fi 2 e(k) -n 2 [6® 6](k) + n\9® 6](k) + ■ 



(4) 



including second order convolutions in 6(k) and S g (k), while ne- 
glecting third and higher order terms. 



2.1 The linear regime 

If 9 and S g are small, then we can drop the second and higher order 
terms from Eq. {4j, and 



s s Jk) = s g (k)-n 2 e(k) 



(5) 



Often it is further assumed that the velocity field comes from 
linear perturbation theory. Then 

9(k) = -/<W*) , (6) 

where / = dlnD/dlna a Si° m 6 dPeebleslll980l) . 

For a population of galaxies, which we denote with a subscript 
g, the linear redshift-space power spectrum can be written 



P s g (k, M ,lm) = <|<S*(M)I 2 >, 

= P gg (k)-2^ 2 P gfl (k)+^P ee (k) 



(7) 
(8) 



where P gg (k) = <|<5 g (k)| 2 >, P g6 (k) = <<5 s (k)0(k)>, P m {k) = <|0(k)| 2 >, 
are the galaxy-galaxy, galaxy-0 and 0-6 power spectra respectively 
for modes k. A subscript shows that a variable is determined from 
the velocity field of the galaxies, which we have assumed is irrota- 
tional and small compared with the real-space distance to the galax- 
ies (see Fig. QJ. In the following we often drop explicitly showing 
the k dependence of these power spectra, for convenience. 



2 MODELLING GALAXY CLUSTERING IN 
REDSHIFT-SPACE 

The redshift-space position of a galaxy differs from its real-space 
position due to its peculiar velocity, 

s = x - j(-(x)z , (1) 
1 The fits are quite sensitive to the satellite fraction, for example. 



2.2 Galaxy bias 

It has long been known that galaxies do not trace the mass, a phe- 
nomenon known as 'bias'. There are good theoretical reasons to 
believe that on large scales th e shape of the power spectra o f galax- 
ies and mass are similar (e.g. |Peebles|[l980l : |Peacoc3ll999h . which 
can be phrased as the assumption of scale-in dependent bias. Un- 
der this approximation. IPekel & Lahavl d 19991) introduced the bias 
relation, 
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Figure 1. Mass density and velocity power spectra recovered from three 
different simulations. Solid lines correspond to the primary simulation used 
in this paper, and for this simulation we show power spectra of the over- 
density (solid squares), velocity divergence (solid circles) and the curl of 
the velocity (solid triangles). For comparison, the dashed line is the corre- 
sponding power spectra from a simulation with the same initial conditions 
and cosmology, run with a PM code (so lower force resolution), and the 
dotted lines corresponds to a simulation with a box half as big (so twice the 
force resolution, 8 times the mass resolution but 1/8 the volume). 



depend on the mass overdensity not the galaxy overdensity. To see 
this, we now consider the normalisation of the power spectrum. 

2.3 Power spectrum normalisation 

The normalisation of power spectra is usually measured using the 
rms fluctuation amplitude, in spheres of radius 8/T 1 Mpc, which 
would have arisen if linear evolution had proceeded to the time of 
observation. For the anisotropic redshift-space power spectrum we 
can define a slightly non-standard, anisotropic normalisation 

°»*M = J ^Wl(k)k 2 Pl(k,rf , (13) 

where Wg(k) is the Fourier transform of the top-hat window func- 
tion of width 8 h Mpc - . Here we have decomposed the standard 
3D integral into a ID integral assuming statistical isotropy of the 
density field, and have performed this integral along particular di- 
rections to the line-of-sight. For each /i, this is equivalent to the 
standard definition of erg for an isotropic power with the same am- 
plitude as P'(k,fi). Substituting Eq. ([§) into this expression, and 
using the fact that P gg , P gB & P se are isotropic, we have that 

°i, g W = (*o"8.m ass ) 2 +2^ 2 (forcr 8-mass )(/o- 8 , mass )+^ 4 (/cr 8 mass ) 2 .(14) 

Because of the \i dependence, large-scale peculiar velocity obser- 
vations provide a measurement of /cr 8 raass that is independent of 
a local linear bias. From observations of the large-scale amplitude 
of M(k,(i) as a function of //, we can hope to measure /cr 8 mass , 
fccr 8 mass , and frrcr 8imass . Eq. dl4t shows in compact form that we 
cannot break these parameter combinations and constrain /, b or 
°"8,mass independently solely using P g (k,/j.), although we can mea- 
sure r. 



<<S(x)<S(x)> (S(x)S g (x)) 
(6 g (x)6(x)) (S g (x)S g (x)} 



<l<Wx)r> 



1 br 

br b 2 



(9) 



where b is the bias factor, and r is the dimensionless correlation 
coefficient between the distributions of mass and galaxies. If r = 1, 
we have a fully deterministic local linear bias relation, 



6 g (x) = M mass (x) 



(10) 



If biasing is not a purely Poisson process, but has a random stochas- 
tic element, then we could h ave r < 1. Measureme nts of r, such as 
summarized in Figure 1 1 o f lSwanson etal] < |2008|) . indicate that r 
is close to 1 on large scales. 

If b = 1, or <5(x) and bS(x) are small, then galaxy velocities 
randomly sample those of the matter in the Universe and we can 
write, 



Pm(k) 
P g e(k) 



fP«m(k) 
-bfrP m , si (k), 



(11) 



where P mass (fe) = (|<5(fc)| 2 ). This is not true for some models of 
galaxy bias, such as the peaks model, which predicts a velocity- 
bias where t he distribution of galaxy velocities does no t match that 
of the mass teegos & Szalay|l995l : |percival et al]2008l) . Under the 
assumptions that lead to Eqns.[TT] the standard way of writing the 
redshift-space galaxy power spectrum is 



P s g (k,fi) = P gg (k)[\+2^/3 + ^Y] 



(12) 



where [} = f lb. By writing the power spectrum in this form we are 
hiding a key feature of linear redshift-space distortions: that they 



2.4 The quasi-linear regime 

In addition to removing higher order terms, the two assumptions 
leading to Eq. ® are that we are dealing with an irrotational ve- 
locity field, and that we are in the distant observer limit (u z <k z) 
both of which we assume to remain true in the quasi-linear regime. 
Such assumptions (among others) are also made in standard Eu- 
lerian perturbation theory djuszkiewicz 198 ll ; Vishniadl 19831; IftvI 
1 19841 ; iGoroff et alJI 19861 ; iMakino et alJI 199a IJain & Bertschingeil 
1994; Bernardeau et al. 2002), Lagrangian perturbation theory 
(Buchert 1989; Hivon et al.lll995l) . renormalised perturbation the- 
ory {Scoccimarro 200T|) and resummed Lagrangian perturbation 
theory (Matsubara 20q3). 

Going beyond the linear assumption leads to additional terms 
in jj 2 , fi 4 & jj 6 in the redshift-space power spectrum compared with 



Eq. <8}. So, 

3 

P^^qlin^^A,^ 2 



(15) 



where A [ + P gB , A 2 t Pgg, and A 3 4- 0. The importanc e of the extra 
higher order terms was emphasised most recently bv lScoccimarrd 
(2004). In addition to these redshift-space effects, the linear theory 
relation between 5 and 6, given in Eq. (6J, will break down in the 
quasi-linear regime, so we should expect the shapes of P gg , P g g and 
Pgg to be different. The relationship between Eqns. ll8l&ll5t can be 
written 



P* g (k,n, qlin) = G(fc,//)^(M. Mn) 

where G(k,ji 2 ) has the property lim^o G(k,/j. 2 ) = 1. 



(16) 
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2.5 The non-linear regime 

The standard model for redshift-space distortions includes a com- 
ponent caused by an uncorrelated velocity dispersion that grows on 
small scales. Such a model is motivated by the idea of "thermal mo- 
tion" of particles in collapsed structures, wh ich causes the Fingers- 
Of-God (FOG) observed in redshift surveys djacksonl 19721) . An ad- 
ditional component comprising uncorrelated particle motions will 
dilute both the galaxy overdensity S g , and the "extra" overdensity 
term caused by the linear distortions, 6. Motivated by numerical 
simulat ions dShefh & Diaferidl200ll ; lHuff et alj|2007l) and the halo 
model dWhitdl200ll ; ISeliakll200lh . we can assume that the centre 
of mass of a halo, around which galaxies orbit, still moves accord- 
ing to (quasi-)linear mot ion. Such a mo del leads to the much-used 
'streaming' models (e.g. Hamilton 1998j) where 



P'(.k,/jM) = F(k,i/)P°(k,ti,ttn) 



(17) 



with F(k,/A 2 ) a function that depends on the distribution of ran- 
dom pair velocities in collapsed objects, which is often written as a 
function of y = kcr, where a is the rms velocity dispersion. In order 
to match behaviour on large scales, we require lim^o F(k,fi 2 ) = 1. 
We have written the equation in this form to highlight the similarity 
with Eq. U6b . Note that this model is constructed by a rather ad- 
hoc splicing of linear, quasi-linear and non-linear behaviour which 
ignores the scale-depen dence of the mapping between rea l and red- 
shift space separations jFishej[T995l : [Scoccimarro 2004], and ref- 
erences therein), while Eq. dl6t was based on the analysis of the 
redshift-space distortions in the quasi-linear limit. 

In general FOG are difficult to model well, and their amplitude 
is strongly dependent on the mean hal o mass and satellite fractio n 
of the population under consideration (White 2001; Seliakj l200ll) . 
Previous work has conce ntrated on models with Gaussian or Expo- 
nential distributions (e.g. ICole et al.ll99l ; |Peacock & Doddsll996l) 
for the pairwise velocity dispersion in configuration space. For an 
Exponential model for the pairwise velocity dispersion in config- 
uration space, we expect a Lorentz damping factor for the power 
spectrum, while the Gaussian dispersion translates to a Gaussian 
damping of the power spectrum 

FExpo„e„,,al(fc,p 2 ) = [ 1 + (fa^) 2 ]" ' , (18) 

^ Gaussian 0c, I?) = eXp [-(lcr/z) 2 ] . (19) 

These terms have the same behaviour to first order. 

The exact form of F(k,/j 2 ), an d the value of cr is strongly 
depen dent on the galaxy population ( Jing & Borner 20041 ; ILTet al.l 
120071) . An alternative approach would be to try to "eliminate" the 
FOG by applying a halo finding algorithm to the sample and manu- 
ally moving galaxies either to halo centres, o r to a spherically sym - 
metric distribution around these centres (e.g. lTegmark e t al. 2004). 
Such approaches tend to mask the fact that the radial "compres- 
sion" is still model dependent, and requires a similar free param- 
eter to cr in Eq. j20\ . This "parameter" controls the probability 
density function for th e distortion of any galaxy in redshift space 
(Reid & Spergel 2008). However, FOG compression does have the 
advantage of including extra information in the analysis from the 
phases, which are used to locate the halos. 

Dealing directly with the FOG is not the same as extending 
the linear model, and P gg , P gB & P m into the non-linear regime, 
as discussed in the previous section. The real-space effect of ran- 
dom thermal motion of galaxies on small scales would lead P gB to 
decrease in amplitude, because of the decoherence of density and 
velocity divergence, while P m increases. We showed in Eq. ( |16t 
that a similar function to F(k,fi 2 ) would be required to include 



quasi-linear behaviour in the redshift-space power spectrum. In this 
case we would not expect that these simple damping models can si- 
multane ously match both the quasi-linea r and fully non-linear be- 
haviour jFisheal 1995l ; IScoccimarro 2004). 

Clearly, if F(k,fi 2 ) has two roles, we must reconsider any im- 
plied or assumed physical meaning in the function. So, for example, 
the physical arguments that the Exponential form sho uld be a bet- 
ter match to the average velocity disp ersion in halos lShetti ll996l ; 
Peacock & Dodds 1996; White 2001) are not applicable. For con- 
sistency with the standard 'streaming model', we continue to refer 
to the function F(k, /j 2 ) as a small-scale velocity dispersion (SSVD) 
model. 



2.6 Combined model 

Following our review of the theory behind redshift-space distor- 
tions presented in Sections [2.1M2.5l we are left with a model of the 
redshift-space galaxy power spectrum 

P s g (k, it) = [P gg (k) + 2v 2 P gB (k) + fi 4 P ee (k)] F(k, S). (20) 

To first order in k, both the Gaussian and Exponential forms for 
F(k,fi 2 ) considered in Section |231 act as an extra yu 6 component 
with amplitude determined by cr. This fitting function can act as an 
approximate correction for both quasi-linear and non-linear terms 
in the mapping between linear real-space and observed redshift- 
space. P sg is independent of galaxy density bias, and is directly 
related to the matter velocity power spectrum, provided there is no 
velocity bias, caused by a mis-match between the distribution of 
galaxy velocities and the distribution of velocities in all matter. 

Eq. d20t could be fit ted to data directly using a Likelihood 
approach ( Hamilton 2000; Tegmark et al. 2004, 2006), resulting in 
observational constraints on P gg , P gB , and P ge . Linking P gg , P gS , 
and Pgg to the matter power spectrum requires further assumptions. 
Following the deterministic, linear local galaxy bias model pre- 
sented in Eq. dlO| >, and the assumption of no velocity bias, the three 
power spectra are related to the matter power spectrum according 
to Eqns.Qj] which leads to a model 



PQcn) = [b 2 + 2f, 2 bf + ^f 2 ]P m .^(k)F(k,n 2 ). 



(21) 



This allows observations to constrain the matter power spectrum 
multiplied by f 2 . We show that this model matches simulation re- 
sults on large-scales in Section|6] If Eq. ( 121b holds, Likelihood fits 
could be used to measure /crg(mass) and icrg(mass) from the mea- 
sured power spectra. Alternatively, we show in the next Section 
that estimators of the matter power spectrum multiplied by f 2 can 
be constructed, based on a Legendre polynomial decomposition of 
the redshift-space power spectrum. 



3 REVISED ESTIMATORS FOR COSMOLOGICAL 
INFORMATION 



We argued in Section l2Tl that large-scale linear redshift-space dis- 
tortions should enable us to measure /tr 8 (mass) independently of 
linear bias. We will now show how we can do so, based on a Leg- 
endre decompositon of the redshift-space power spectrum. P" g (k,{i) 
can be decomposed into Legendre polynomials L e (ji) to give mul- 
tipole moments 



P s e (.k) : 



21 + 1 



diiP s (k,[i)L ( {n) 



(22) 
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The first three even Legendre polynomials are L (jj) = 1, L 2 Qi) = 
(3yu 2 - l)/2 and L 4 (p) = (35yu 4 - 30yi/ 2 + 3)/8. For redshift-space 
distortions along one axis, the expansion of the linear model for 
P g (k,fi) (Eq.[8]l into the first three multipole moments gives 



C p° m > 




■ 1 


2/3 


1/5 




( p sg ( k ) ) 




P s 2 (k) 







4/3 


4/7 




P g e(k) 


(23) 


{ P' 4 (k) j 




I o 





8/35 


) 


{ Peo(k) J 





so we see that the monopole, quadrupole and hexadecapole com- 
pletely characterise P s (k) in the linear limit. 



If we assume that P„ 



> have different forms, 



we can invert Eq. d23b to recover them individually from Pg P- 



then 
, and 



( P ss (k) 1 




f 1 


-1/2 


3/8 




( p o(k) ) 




P gS (k) 







3/4 


-15/8 




P s 2 (k) 


(24) 


, Pw(k) ) 




( o 





35/8 


7 


{ Pl(k) J 





In terms of fi, P gg , P g e, and Pee can be considered as revised mo- 
ments by substituting the expressions for P^, P s 2 , & P s 4 into this 
expression. 

The inclusion of a SSVD model given by Eq. {78} or Eq. d!9t 
in the analysis of this statistic is considered in lCole et alj dl995T) : 
the integrals in Eq. d22t need to be recalculated as in Section [5^21 
This would enable us to use the monopole, quadrupole and hexade- 
capole to extract the same information from simulations as using 
the spherically averaged power spectra described in Section|5] 

Following Eq. d 1 2b . measur ements of P^/P s n have been previ- 
ously used to measure j3 through dCole etal.ll 19941) 



P s 2 (k) 

pm 



i + f/j+is 2 



(25) 



This only involves the "lowest order" moments in \i (P^ and PI) and 
it has been argued that this makes it less sensitive to noise. Expand- 
ing in multipoles can also simplify the definition of simple survey 
limits. The expansion of thi s formulae to inc lude a SSVD model 
has also been considered bv lColeetalJdl995» . The quadrupole to 
monopole ratio is expected to have this limiting behaviour on large- 
scales, and can therefore be used to measure fi. 

We now show that we can use P S Q and PI, to derive an esti- 
mator of the matter power spectrum multiplied by f 2 , which leads 
to cosmological constraints that are independent of a local galaxy 
bias. For a local linear bias, in the absence of velocity bias, P m and 
(P g e) 2 /Pgg both have normalisation / 2 o"g(mass) and shape match- 
ing that of the matter density power spectrum in the linear limit. Re- 
sults from simulations, presented in Section [6] are consistent with 
this claim. W can therefore define a new estimator P, such that 
P = (P g e) 2 /P sg , and P = Pge on large scales. 

In terms of the Legendre decomposition, we can write two 



functions of PI 



P\ and P\, 



which should have the same large-scale 



normalisation that is independent of a local linear bias. In the linear 
regime, from Eq. J24t these functions are 

(26) 



P(k) 



\-P s 



-P s + 

2 2 



J ps 
S r 4 



P(k) 



35 „. 
— Pi 



(27) 



Because we have two equations for P, we can eliminate P s 4 , leaving 
a quadratic equation in P 

24[P(k)] 2 - 35(7^ + P 2 )P{k) + ^r(^) 2 = . (28) 

lo 

Solving this equation gives 



P(k) = -[5(7P^ + Pl)- 



V35[35(^) 2 + WP^Pl 



1(P 2 ) 2 ] 112 ] 



,(29) 



which offers a mechanism for removing a local bias dependence 
from quadrupole and monopole measurements in the distant ob- 
server limit. The result is a power spectrum whose large-scale shape 
should match that of the mass, and normalisation should be f 2 
times that of the mass density power spectrum. 

We can extend this model to the quasi-linear regime by in- 
cluding the model for small-scale velocity dispersion. As we ar- 
gued previously, in the quasi-linear regime, this should be consid- 
ered to be a fitting formula for quasi-linear distortions, and is not 
physically motivated. In the following we only consider a Gaus- 
sian SSVD model, because this is favoured by simulation results 
presented in the next section. A similar formula could be calcu- 
lated for the Exponential SSVD model. In the linear-Gaussian 
model, we can simplify the equations by defining T„(y) = y([n + 
l]/2,y 2 )/(2/ +l )sothat 



P\ 



f(3r 2 -r ) 



2r 2 
15T 4 - 5r 2 



f(3r 6 -r 4 ) 



p, 
P 



(30) 



where we have defined _P X 



^jp~gg 



P for notational convenience. 



This equation reduces to the relevant expression in Eq. (29J in the 
limit as y — > 0. As outlined above for the linear case, these two 
equations can be manipulated to eliminate P gg and calculate an ex- 
pression for P, which is a quadratic in P S Q and P,, as in Eq. l !29t 
and can be easily solved, although it now depends on T„(y) in a 
complicated way. The linear and damped versions of this estima- 
tor are compared with the results from our numerical simulation in 
SectionlMl 



4 SIMULATION 

To investigate redshift-space clustering further we have used a 
large, high-resolution N-body simulation which is well suited to 
probing P g (k,/i) in the quasi-linear regime. The cosmology was 
of the ACDM family with Cl m = 0.25 = I - Cl A , Q b = 0.043, 
h = 0.72, n s = 0.97 and erg = 0.8. The linear theory power spec- 
trum for the initial conditions was computed by evolution of the 
coupled E instein, fluid an d Boltzmann equations usin g the code de- 
scribed in lWhite &s"cottl jl995b . ISeliak et alj d2003l) find tha t this 
code agrees well with CMBfast ( Seljak & Zaldarriaga 1996). The 
simulation employed 1024 3 particles of mass 10" /T 1 M B in a pe- 
riodic cube of side 1 h~'Gpc using a TreePM code dWhitell2002l) 
with a Plummer-equivalent softening length of 35 h~' kpc (comov- 
ing). A detailed comparison of this Tr eeP M code with other co des 
can be found in lHeitmann et alj d2008h and lEvrard et al 1 d2008h . 

In order to test the convergence of this simulation, we compare 
present-day density and velocity power spectra against alternative 
simulations in Fig. [JJ We compare with a simulation run from the 
same initial conditions, but using a particle-mesh code which has 
lower force resolution. We also compare against a simulation with 
half the box size, i.e. double the force resolution over 1/8 of the 
volume. The density power spectra recovered from the simulations 
do not reveal any significant trends with force resolution. There is 
weak evidence for smaller Vxu and velocity divergence with higher 
force resolution, but this is not significant for our purposes. 

The offset between density and velocity divergence power 
spectra seen in Fig.[JJon large scales is caused by the factor f 2 as 
described in Section |2TI The difference on smaller scales is anal- 
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ysed in detail later. Our assumption that the velocity field is curl- 
free is validated by these simulations, which show that |V x u\ k <K 
|V • u| t for the scales of interest: k < 0.2 h Mpc~' . 

From the z = output we gene rate a halo catal og using the 
Friends-of-Friends (FoF) algorithm (Davis et"ail 19851) with a link- 
ing length of 0.168 times the mean inter-particle spacing. This pro- 
cedure partitions the particles into equivalence classes, by linking 
together all particle pairs separated by less than a distance b. The 
halos correspond roughly to particles with p > 3/(27rfe 3 ) =; 100 
times the background density. We will present results calculated 
from all particles and from catalogues containing all particles in 
halos more massive than log 10 M = 12.5, 13 and 13.5 where M is 
in units of h~ l M Q . The small-scale velocity dispersions of the parti- 
cles will match the distribution of mass in each halo. We do not in- 
clude further small-scale effects caused by the galaxies inside each 
halo not Poisson sampling the mass. Our intention in this paper is 
not to try to model the non-linear regime perfectly, but to test the 
effect of halo mass selection on the linear and quasi-linear regimes. 
We therefore adopt this simple halo occupation distribution in or- 
der to reduce the shot noise: by using all particles in each halo we 
optimise our determination of the (quasi-)linear halo velocity. 

We cannot easily replicate Fig. Q] for our halo-particle cata- 
logues because of the difficulty of interpolating the velocity field 
between the sparsely sampled galaxies. Instead, in the rest of 
this paper, we directly analyse redshift-space distortions through 
P g (k,p). We will consider two ways of analysing P g (k,p), either 
using spherically averaged power spectra with different redshift- 
space distortion components, or by decomposing P g (k,p) into Leg- 
endre polynomials. In the simulations we determine PI in each k- 
bin both by directly integrating against p as in Eq. \22i and by a 
least-squares fit of P g (k,p) to a sum of Legendre polynomials. The 
results are consistent. Before presenting the results from simula- 
tions, we first consider the theory behind decomposing simulation 
results into spherically averaged power spectra. 



5 SPHERICALLY AVERAGED POWER SPECTRA FROM 
SIMULATIONS 

In this section we introduce the concept of spherically averaged 
power spectra calculated from simulations where we propagate 
redshift-space distortions along multiple axes. Clearly, such power 
spectra will not match observational results, where redshift-space 
distortions are only radial. However, they do provide a convenient 
and simple way to explore the components of Eq. \20\ using sim- 
ulations. For a sample of galaxies where we have applied redshift- 
space distortions along multiple axes, for each k, Eq. j20\ will still 
hold for some value of p. The spherically averaged power spectra 
act as /^-dependent moments of P g (k,p), and this method does not 
alter the physics that we are testing. We are simply summing these 
modes over different combinations of p. When analysing simula- 
tions, we apply the plane-parallel decomposition to both the simu- 
lation data and the analysis method. Consequently, this should not 
be considered an approximation as we do not relate the results to 
actual surveys, but only use simulations to investigate the physics. 



For two axes with redshift-space distortions, we can take v to be the 
cosine of the angle to the direction where there are no distortions 
and integrate so that, in the linear limit, 



f 

Jo 



dv ([5M + 0(k) 



(32) 



For three axes, whatever direction we take our fc-vector, we see 
redshift-space distortion of the overdensity is as in Eq. (5} with 
p = 1. Solving the integrals in Eqns.OTT) & J32b . and substi- 
tuting in the relations P gg (k) = (\S g (k)\ 2 ), P g9 (k) = (S g (k)6(k)) 
& Pgg(k) = (\6(k)\ 2 ), leads to the following spherically averaged 
power spectra, 
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(33) 



Here P^ axes (k) includes redshift-space distortions along ( axes. Us- 
ing any three of P, ! axes (fc), we can reconstruct P gg , P g9 , and P m by 
solving the corresponding linear equations. 

5.2 Including small-scale velocity dispersion 

If we include a SSVD damping model then, for a power spectrum 
with redshift-space distortions along one axis, 



,(*)= ( dp ([6 g (k)+p 2 G(k)f)F(k,p 2 ) 



(34) 



For redshift-space distortions along two axes, we take v to be the 
cosine of the angle to the direction where there are no distortions, 
as in Eq. d32b . We can consider that the power spectrum is damped 
by the factor F(k, 1 - v 2 ) so that, 



Jk) 



f dv ([S g (k) + (1 - V 2 )6(k)?F(k, 1 - v 2 )} 
*Jo 



(35) 



For redshift-space distortions along three axes, whatever direction 
k/k we choose, we see redshift-space distortions as if we were 
looking along the line of sight. So 



Jk) = ([SJk) + 8(k)YF(k,l)) 



(36) 



If we assume an Gaussian model for the SSVD damping, then 
Eq. d34t can be solved using the y function, 



y(a, x) '■ 



f 

Jo 



dte-'f 



(37) 



and the factors 1, 2/3 and 1/5 in the expansion of PJ . (k) in 
Eq. (33} become y(l/2, y 2 )/(2y), y(3/2,y 2 )/y 3 and y(5/2,f)/(2y s ) 
respectively. The original factors are recovered in the limit y — > 
0. Alternative expressions can be derive d using error functions 
(Cole et alJll995l : |Peacock & Dodds|[l996h . Eq. (35} can be solved 
similarly using the imaginary error function. 

If instead we assume an Exponential model for the small-scale 
velocity dispersion, then the inte gral in Eq (341 ) can be analytically 
solved using the tan -1 function (Cole et aljl995h . while the integral 
in Eq. (35} is trivial. 



5.1 Linear distortions 

For redshift-space distortions along one axis, in the linear regime 
where Eq. (5} holds, we have 



W*)= f d v {[s g (k)+p 2 e(k)f) 



(3D 



6 SIMULATION RESULTS 

6.1 Fits to the spherically averaged power 

The first question we wish to address is "How well does the linear 
model with or without SSVD model (Eg. 121b recover the monopole 
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Figure 2. Recovered power spectra from our simulation, with redshift-space distortions along one axis of the simulation box (solid circles). Power spectra 
are shown for all mass (lower data), and for the three halo catalogues described in the text (upper data, with the more biased power spectra corresponding to 
larger mass thresholds). These data were fitted for k < 0.05 h Mpc~' using the model given by Eq. 12 1 i . Ratios between model and data are shown in the lower 
panels, and panels from left to right are for different SSVD models. See Section[5]2]for details. 



power P s = P s lilxis ?' This comparison is shown in Fig. [2] The two 
possible free parameters, f/b and cr, were fitted on very large scales 
k < 0.05 ZzMpc -1 . The amplitude of the power spectrum on these 
scales is close to the expected value for all data, following linear 
theory (Section r2.lt and given the simulation value of /?. Where no 
SSVD model is present, we see that the shape of the model power 
spectrum does not match the data, even on these large scales. In- 
cluding either the Gaussian or Exponential SSVD model allows a 
better match of the shape. The goodness-of-fit is approximately the 
same for both of these models, which is not surprising as they have 
the same asymptotic behaviour in the limit y — > 0. High values 
of cr ~ 500-600 kms~' are required to enable this match, which 
fails for k > 0. 1 h Mpc" 1 . While the SSVD model was introduced to 
correct the high-fc behaviour of redshift-space power spectra, it is 
clearly having an effect on scales 0.05 < k < 0.1 ZzMpc -1 , where it 
is attempting to correct both the quasi-linear redshift-space effects 
and the deterministic linear bias assumption. 

The second question that we wish to answer is "How well 
does the model of Eq. ((20} where P gg , P gB & P m are arbitrary, with 
or without a SSVD model, fit the simulation results?". The power 
spectra, P\ . are shown by the solid circles in the upper panels of 
Fig. [3] On large-scales, k < 0.1 /zMpc~', the data are well fitted by 
the model of Eq. < |20| > without the inclusion of a SSVD component. 
As we move to scales k > 0. 1 h Mpc - ', we need to include a SSVD 
model in order to match the differences between Pi . , Pi , and 

! axis' 2 axes' 

faxes' wn i cn result from the ji dependence of P g (k,/j). Ratios be- 
tween model and data for Pi , PI , and Pi are shown in the 

1 axis ' 2 axes ' 3 axes 

lower panels of this plot. Although the SSVD models significantly 
help with the fits to the observed power spectra, it is clear that the 
shapes cannot be matched and that a compromise is being reached 
in the fits, where the model crosses the data at k ~ 0.6 ft Mpc -1 . 
The best-fit values of cr vary with the range of scales fitted confirm- 
ing this picture. On these quasi-linear scales, the Gaussian SSVD 



model allows a slightly better fit to the data, although it is not clear 
that there is a physical motivation behind this. The Exponential 
model is a better fit on smaller scales, whi ch is physically moti- 
vated jPeacock & D odds 1996; White 200lb . 



6.2 Extracting density and velocity power spectra 

The third question that we want to use the simulations to answer 
is "How are the power spectra P gg , P gf) & P m related to each other 
and to the matter power spectrum?". To fit for P gg , P gS & P SB we 
use spherically averaged power spectra calculated assuming that 
we have redshift-space distortions along 0, 1, 2 or 3 axes. Without 
including a SSVD model, we could use any three of Pi , Pi . , 

t> ' J axes ' 1 axis ' 

faxes' and ^3axes t0 S ° lve Ec l- f° r P gg> P S» 3Ild P M' In 0rder 

to simplify the fit, we assumed P„„ = Pi , and used Pi . and 

r J » iS axes ' 1 axis 

^2axes t0 s °l ve f° r Pge & Pm- For fits including a SSVD model, we 
used all four power spectra to solve for cr, P gg , P ge and P m . 

The power spectra, P gg , P g g/f, Pgg/f 2 , recovered from the fits 
to P?„j- are plotted in the upper panels in Fig. [4] The power spectra 
divided by Pmass are shown in the lower panels to highlight varia- 
tions in shape of these different power spectra. The residuals be- 
tween model and simulation data for Pi^ were shown in Fig. [3] 
On the largest scales k < 0.05 h Mpc~' , the data are consistent with 
the hypothesis that P m ^ f 2 P lwdss , and P 2 gg - f 2 P gg P mdsi . Where we 
do not include a SSVD model, there is weak evidence that the am- 
plitude of Pg e I 'f 2 is lower than the amplitude of the matter power 
spectrum on scales k < 0.05 /jMpc -1 , consistent with a 10% veloc- 
ity bias, although it should be noted that the data are noisy here. 

It is clear that the shapes of the velocity and overdensity power 
spectra do not match on scales k > 0.05 /iMpc~', even when ei- 
ther Gaussian or Exponential SSVD models are included. Devia- 
tions betwe en the shapes of the se power spectra were also noted in 
Figure 6 of lScoccimarrol d2004l) . and were physically motivated by 
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k / h Mpc 1 k / h Mpc 1 k / h Mpc 1 

0.01 0.02 0.05 0.1 0.2 0.5 0.01 0.02 0.05 0.1 0.2 0.5 0.01 0.02 0.05 0.1 0.2 0.5 




log 1D k / h Mpc 1 l°£io k / h Mpc 1 l°Sio k / h Mpc 



Figure 3. P j . (k) calculated from the simulation are shown by the solid circles - these are the same data in the upper panels of Fig. [2] The lowest amplitude 
power spectrum corresponds to the mass, while the other three are from halo catalogues, as described in the text, where the amplitude is an increasing function 
of mass limit. For comparison we also plot the model given by Eqns.[34] where the free parameters fitted are P gg , P g g, Pgg, and <x where a SSVD model is 
included (solid lines). To determine these models, we have fitted axis , ^? ax(:s > an d f° r k < 0.7 /i Mpc , where the data were equally weighted in log k. 

In order to highlight differences, we plot the model for Pj • divided by the data in the lower panels (solid lines). Similar lines are also plotted for Pi 
(dashed), andP^ axes (dotted). 



k / h Mpc 1 k / h Mpc 1 k / h Mpc 



0.01 0.02 0.05 0.1 0.2 0.5 0.01 0.02 0.05 0.1 0.2 0.5 0.01 0.02 0.05 0.1 0.2 0.5 
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Figure 4. Upper panels: P gg (dotted lines), P g gl f (dashed lines), Pgg/f 2 (solid lines) recovered by fitting Pg axes , Pi^i faxes' ^3 axes ^ or * < 0J hMpc^ 1 . 
Four power spectra are plotted, corresponding to the mass and the halo catalogues. The amplitude of P gg and P g g is lowest for the mass, and is an increasing 
function of the halo mass limit. For Pgg, the deviation from the shape of P gg is weakest for the mass, and is an increasing function of halo mass limit. The 
lower panels show these power spectra divided by P gg for the mass to highlight deviations between the shapes of the power spectra. 
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Figure 5. For P g g, we can remove a deterministic linear galaxy density 
bias dependence by dividing by ^JP gg - If the galaxy density bias is well 
described by a local linear model and there is no velocity bias, then we 
expect P 2 g0 /P gg — ^mass* so we divide P g g by a further factor of VP mass 
to highlight deviations from this model. The solid line was calculated from 
a fit with no small-scale dispersion correction, the dotted line assumed an 
Exponential model, and the dashed line a Gaussian model, both with single 
scale-independent variance. 



the arguments in Section |2l4l The turn-off from linear behaviour is 
stronger for P se than P gg for any of the catalogues, which is con- 
sistent with the quasi-linear particle velocities leading the displace- 
ments, which is an integral over the velocities. 

As we discuss in Section [3] most of the cosmological signal 
from observations comes from the cross correlation between over- 
density and velocity fields. So, the forth question that we wish to 
use simulations to answer is "Is P 2 g /(P gg P matis ) = 1 a valid model 
on large scales?". To answer this, we plot the power spectrum 
combination P 2 g /(P gg P m ass) in Fig. [5] On large scales, this statis- 
tic — > 1, validating our hypothesis. Comparison with Fig. [4] sug- 
gests that it is less biased than P m as a tracer of / 2 <Tg(mass) on 
these scales, as the deviations from P 2 g /(P gg P m ass) = 1 are less 
than 10%. The large-scale "plateau" can be extended by assuming 
a Gaussian model for the small-scale velocity dispersion, provided 
the fit is not extended to large k. For the plot, we only fitted to 
k < 0.4/iMpc -1 . The recovered values of cr for Gaussian model 
of scale-independent velocity dispersion correction are 310kms~', 
360kms~', 390 kms~', 470 kms~', for the mass, and the three halo 
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Figure 6. As Fig. [4] but now allowing the variance of the assumed SSVD 

model, cr to vary with scale. The resulting models of P\ , Pi , and 

J ° 1 axis' 2 axes ' 

Pj are an almost perfect fit to the data. This should be expected as we 
are fitting 4 data points with a model with 4 parameters at each lvalue. 



catalogues. For the Exponential model, the corresponding numbers 
are 280 kms" 1 , 300 kms" 1 , 340 kms" 1 , 410 kms" 1 . If we extend the 
fit to smaller scales, then cr increases for both models and the re- 
sults from Exponential and Gaussian models become more similar. 
In addition, the plateau is reduced in size. 

The sharp decrease in P g!) seen at k ~ 0. 1 h Mpc~ 1 is consistent 
with a model that includes both the effect of FOGs and the infall 
model. The inflection marks the point where the highest galaxies 
velocities change from being associated with the largest scales (due 
to infall), to instead be linked to the smallest scales (due to FOG). 
The quasi-linear nature of infall on small scales will also contribute 
to this transition. The scale at which the inflection occurs is affected 
by the SSVD model included. The Gaussian SSVD model acts as a 
slightly more accurate fitting function compared with the exponen- 
tial model, extending the plateau of where the overdensity-velocity 
cross power spectrum matches the shape of P gg to k < 0. 1 h Mpc -1 . 
This might reflect the contribution of quasi-linear effects to the 
FOG, where the Exponential model was physically motivated. We 
now investigate the SSVD model further. 



6.3 Investigating SSVD fits 

The fifth question that we wish to address using simulations is 
"How well do the Gaussian and Exponential SSVD models in 
Eq. d201 > extend the fit to small scales?". This is clearly a diffi- 
cult question to answer, but we have tried to quantify the effect, 
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and determine the scales where the SSVD model is important, by 
allowing the free parameter, a in the Gaussian and Exponential 
SSVD models to vary as a function of scale. The power spectra 
P gg , P g e/f & Peel 'f 2 , calculated assuming the Exponential SSVD 
model are shown in Fig. [6] which has the same format as Fig. [4] 
Allowing cr to vary allows a virtually exact fit to P^, faxes' an ^ 
P 3axes- For k < al h M P C ~' - Fig.[3]shows that no SSVD model is re- 
quired to fit the observed power spectra, and there is no constraint 
on cr. For k > 0. 1 AMpc -1 the best fit value of cr rises from at 
k = 0.1 hMpc~ l to 295 kms~' at k = 1 ftMpc" 1 for the matter, and 
to 330kms -1 , 390 kms -1 , and 470 kms~' for the halo catalogues 
with mass > 10 1Z5 , 10 13 & 10 13S h~' M Q respectively. A similar 
trend is seen for the Gaussian model, although there is more noise, 
perhaps indicating that the parameters P gg , P gS , & Pgg and cr fit- 
ted for each k are more degenerate in this model. This would also 
explain why the Gaussian model with a single value of cr for all 
scales can provide a slightly better fit to the data shown in Fig. [3] 
the solution can move along the Likelihood degeneracy. 



6.4 Testing our new estimator 

The final question that we wish to answer from simulations is "How 
well do the power spectra constructed from our new estimators 
given in Eqns. i24\ & l |30t , recover f 2 P milss (k)7"- The power spec- 
tra are compared in Fig. [7] Even with only P^ and PL we see that 
we can recover an unbiased estimate of / 2 _P miiss (£) on large-scales. 
The inclusion of a Gaussian SSVD model is able to correctly model 
the quasi-linear behaviour of this estimator for k < 0.2 h Mpc~' for 
the mass, a limit that decreases to k < 0.1/zMpc~' for the most 
massive halo catalogue. The value of cr can therefore be treated 
as a free "nuisance" parameter, and could be fitted to the estima- 
tor. In order to demonstrate consistency we actually plot the esti- 
mator as in Eq. l |30| l, with cr fitted to P' . , Pi„„, and Pi for 

1 1 axis I. axes j axes 

k < 0.4/zMpc~'. In practice, for a given linear matter power spec- 
trum, one could calculate the estimator for a series of values of cr 
and marginalise over this parameter when fitting to the model. 

Comparing the estimator of Eq. J30| > to the non-linear matter 
power spectrum would result in a better fit to the data, because the 
behaviour of P and the non-linear matter power spectrum are simi- 
lar. We consider that this would be artificial, because we do not ex- 
pect such a match from the theory presented in Section|2] With the 
Gaussian damping term, our new estimator is clearly a complicated 
function of P^ and PI, and it would be difficult to accurately prop- 
agate standard analytic error estimates for power spectrum modes 
to this function. Given the precision to which future surveys would 
be able to measure this statistic, it would be better to base error 
analyses on the results from mock catalogues. 
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7 DISCUSSION 

Redshift-space distortions encode key information about the build- 
up of cosmological structure. In the linear regime, the theory be- 
hind the den sity enhancement was developed over 20 years ago 
jKaiserlll987l) . Following this development, most observational 
studies have focused on measuring f}. One method to do this is to 
measure the quadrupole-to-monopole ratio, which depends on [} in 
a known way (Eq. 125b. By reviewing the theory of redshift-space 
distortions in Section|2] we argued that we should be able to obtain 
a bias independent cosmological constraint on /<r 8 (mass) from the 
large-scale redshift-space distortions. This statistic, proportional to 
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Figure 7. The power spectrum given by Eq. i29\ divided by the real-space 
linear power spectrum multiplied by f 2 (open triangles). On large scales 
k < 0.05 h Mpc~' this tends to 1 as expected. For comparison, the solid 
triangles give the value of the estimator given by Eq. )30t . where we include 
a Gaussian component equivalent to a SSVD model, with variance fitted to 
P s , . , PI , and P s , (see text for details). This demonstrates that a 
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Gaussian SSVD model on scales k < 0.4ftMpc~' can be treated as a fitting 
formulae for the quasi-linear turn-off from the expected linear behaviour. 



dD/d\na with a proportionality constant depending on the ampli- 
tude of fluct uations at early times, provides an excellent test of DE 
models (e.g. lSong & PercivalH2008l) . 

In Section [2~4l & l2.5l we set out to extend the theory into the 
quasi- and non-linear regimes, to investigate the smallest scales on 
which we can easily recover cosmological information using this 
physical process. In the linear regime, the power spectrum can be 
decomposed into terms dependent on fi°, fi 2 and /u 4 . As we push 
to quasi-linear scales we should expect this behaviour to become 
more complicated with extra jj 6 and higher order terms becoming 
important. We should also expect the simple relations bP gB = fP gg 
and b 2 P e g = f 2 P gg to break down. Previous work has introduced 
a "streaming" model for non-linear redshift-space distortions, in- 
cluding a damping term spliced together with linear theory. In the 
limit k — > 0, the standard Gaussian or Exponential forms given in 
Eqns. i l 1 8t & l !19t for this damping function both reduce to provid- 
ing an extra fi 6 term with amplitude dependent on cr. We might 
therefore expect these terms can also be used to match the ex- 
pected quasi-linear n dependence of the power spectrum. In this 
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case, the same function might not match all quasi-linear and non- 
linear scales. 

Using a large numerical simulation we have set out to test this 
theory. To facilitate this, we have introduced a novel approach to the 
analysis of the redshift-space power spectrum measured from sim- 
ulations by allowing redshift-space distortions to be included along 
multiple axes. Power spectra calculated from density fields includ- 
ing redshift-space distortions along multiple axes act as moments 
of P s g (k,fi), allowing us to decompose into the expected pr and 
fi 4 components. On large scales we see no evidence for strong ve- 
locity bias, and the velocity-velocity power spectrum has an ampli- 
tude f 2 times that of the mass to within 10% (right panel of Fig.[4]l. 
The density-velocity power spectrum, from which most informa- 
tion is obtained from surveys is less affected by this velocity bias. 
Fig. [2] shows that the relations bP g0 = fP gg and b 2 P m = f 2 P gg 
break down on scales k > 0. 1 AMpc -1 . The right panel of Fig. [3] 
shows that // and higher order terms need to be included for 
k > 0.2/iMpc~'. The Exponential and Gaussian models enable a 
better fit for 0.2 < k < OAh Mpc~', but neither is perfect on all of 
these scales, although the Gaussian model does slightly better than 
the Exponential model in the cases we illustrate. Alternative mod- 
els which reduce to yu 6 behaviour in the limit k — > might perform 
better, but this is not necessary for our purposes: we have shown 
that the Gaussian model extends the range of the fit. 

Based on this analysis, we have devised a new estimator of 
the matter power spectrum P using the monopole and quadrupole 
power spectra. This matches the large-scale shape of the matter 
power spectrum with amplitude multiplied by f 2 , so we can mea- 
sure /crg(mass). Including the Gaussian damping model allows us 
to extend the match with the simulation data for k<0.2hMpc~ l . 
The match between the shape of the expected linear matter power 
spectrum and P should allow a test of the scales over which we can 
measure /cr 8 (mass). Rather than use such an estimator, it would 
also be possible to perform a maximum likelihood fit of Eq. d21t to 
the measured redshift-space power spectrum, with /crgCmass) and 
£>cr 8 (mass) as the free parameters. 

Our analysis has assumed the distant observer limit, a stan- 
dard approximation used in many analyses. It would be straight- 
forward to extend this work to methods that incorporate a proper 
radial-angular split allowing for standard survey geometries. An 
additional extension would be to perform such an analysis in con- 
figuration space: again, it would also be possible the extend the 
ideas behind our estimator to analyse the correlation function. 
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